On the single mode approximation in spinor-1 atomic condensate 
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We investigate the validity conditions of the single mode approximation (SMA) in spinor-1 atomic 
condensate when effects due to residual magnetic fields are negligible. For atomic interactions of 
the ferromagnetic type, the SMA is shown to be exact, with a mode function different from what is 
commonly used. However, the quantitative deviation is small under current experimental conditions 
(for *^Rb atoms). For anti-ferromagnetic interactions, we find that the SMA becomes invalid in 
general. The differences among the mean field mode functions for the three spin components are 
shown to depend strongly on the system magnetization. Our results can be important for studies of 
beyond mean field quantum correlations, such as fragmentation, spin squeezing, and multi-partite 
entanglement. 
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Trapped atomic quantum gases have provided a re- 
markable testing ground for quantum many-body theory 
101 . Since the discovery of the first atomic Bose- Einstein 
condensate ||[, mean field theory has been apphed with 
great success to these systems, arguably because 1) low 
energy atom-atom interaction can be simply parameter- 
ized by a s-wave scattering length age with atoms behave 
as hard spheres of effective radii asc] and 2) most cur- 
rent atomic gases are dilute with densities n satisfying 
nOgc <C f 1^. Increasingly, theoretical and experimental 
attentions are directed towards beyond mean field effects. 
In this regard, spinor-1 atomic condensates have become 
a proto-type system for many recent studies ||, ||, 
Several interesting results have already been obtained, 
e.g. multi-particle and continuous variable type entan- 
glement Q , spin- mixing j^] , spinor four- wave mixing , 
and super and coherent fragmentation [Q. The single 
mode approximation (SMA) is often adopted for these 
studies when a mean-field approach with a vectorial or- 
der parameter becomes inappropriate |jl^, |l^, 12 1. Be- 
yond mean field quantum effects have been found both 
when there is no external fields 0]and when there is an 
external magnetic or optical field |13| , p^ , [Tsl , |l6|, |l^ . To 
justify the use of the SMA, earlier studies often compared 
with solutions of the coupled Gross-Pitaevskii (GP) equa- 
tion for the different spin components and enforced an 
upper limit on the number of atoms |l^ . While there 
is not a generally adopted limit, it is typically estimated 
that N should be less than 10"*, a rather small number 
for current experiments. 

In this paper, we investigate the validity conditions of 
the SMA in spinor-1 atom condensate ||, ||. Our initial 
aim was to provide a reliable thermodynamic phase dia- 
gram for a trapped spinor-1 atomic gas | |l^ . Surprisingly, 
interesting zero temperature results from the coupled GP 
equations reveal intricate relationships of the mode func- 
tions for the three spin components due to the constraint 
on the system magnetization. 

We consider a spinor-1 atomic condensate in the ab- 
sence of an external magnetic field. As partitioned by 
Law et al. pi, the system Hamiltonian, H, separates 



into a symmetric part (under spin exchange) 



drhlC^p^p+'^^i^l^p^!^, (1) 



with C 



Q/3 



7i^V^/2Af -t- V^xt, and an asymmetric part 



(2) 



where '^a (a — 0,±) denotes the annihilation field op- 
erator for the a-th component. F^^^^ y z are the spin 1 
matrix representation, and a summation over repeated 
indices is assumed in Eqs. (|^) and (|^). The external 
trapping potential V^xt(?^) is spin-independent as in an 
far off-resonant optical dipole force trap (FORT) which 
makes atomic spinor degrees of freedom completely acces- 
sible. The pair interaction coefficients are cq = 47r7i^(ao-l- 
2a2)/3A/ and C2 = 47r?i^(a2 — ao)/3M, with ao (02) the 
s-wave scattering length for two spin-1 atoms in the com- 
bined symmetric channel of total spin (2). The only 
state changing collision in Eq. (H) occurs through the 
coupling \E'Q5'g^'+5'_ + h.c, which conserves the system 
magnetization = J df{Fz) / df[^'\.\E'+ - ^E'^l*.]. 
Al-changing inelastic ( "bad" ) collisions occur at a much 
longer time scale as compared with a condensate's typ- 
ical lifetime, therefore are excluded here as in all pre- 
vious studies. Although the real time dynamics gov- 
erned by Hs + Ha conserves the total atom number 
N ^ J dr[vE'^*+ + ^I'^o + ^l^*-] and M, the ground 
state obtained from a global minimization of Hs + Ha 
is not automatically guaranteed to have the same N and 
M. We therefore introduce separate Lagrange multipli- 
ers B to guarantee the conservation of Ai and the chem- 
ical potential /x to conserve N. The ground state is then 
determined by a minimization of the free energy func- 
tional — Hs + Ha — fJ-N — BM. . Mathematically, this 
task turns out to be highly nontrivial. In fact, most previ- 
ous discussions on spinor-1 condensates did not minimize 
H under the constraint of a conserved M. Therefore, 
their resulting ground states are the global ground states 
that can only be reached if the system can coherently ad- 
just its initial M. value. Such a situation is inconsistent 
with current experiments. 
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FIG. 1: The Ml dependence of A'o in the ferromagnetic case. 
The solid line shows No/N = (1 - M^/N')/2, while the dots 
are numerical results. The agreement is remarkable. 

One of the strongest physics support for the SMA 
comes from the fact that oq ^ 02 for a spinor-1 (^^Rb) 
condensate. This gives rise to |c2| ^ |co| [|[ Thus 
Ha is much smaller as compared with Hs, and can be 
considered as a perturbation by assuming the SMA 



a = 0,±, 



(3) 



i.e. with a common mode function (j^sMAir") (normal- 
ized to 1). The Fock state boson operators satisfy 

[o-a, a. 



7J 



] — 0. (/>sma(?^) is determined from 



Hs alone (without Ha) according to 



2M 



coN\(f>. 



SMA I 



It shares similar physics of the often used spin-charge sep- 
aration in condensate matter systems. Since its introduc- 
tion, the SMA has been used frequently ||, ||, |§7 0, p^ . 
Notable exception is the work by Ueda [20|, who went 
beyond the SMA by studying a translational invariant 
system with the use of a plane wave basis. Correlations 
between spatial and spinor degrees of freedom were then 
shown to lead to effects associated with density waves 
and spin waves. For a trapped system as studied here, 
the use of a plane wave basis becomes inappropriate. 
The same SMA is sometimes also used in a spin 1/2 

, I2I. This 



system by assuming 0o(^ = '/'i(^) 1 21 



is less critical as the resulting Hamiltonian oc J|reniains 
of the same symmetry group in the Schwinger boson rep- 
resentation, although with a different coefhcient and the 
presence of additional linear terms in J^. The validity of 
the SMA in this case has been tested recently using the 
rigorous positive P-approach p^ . 

For a spinor-1 condensate, however, complications 
arise when spin component mode functions are taken to 
be different. The effective Hamiltonian thus obtained 
contains no angular momentum symmetry at all in its 
corresponding Schwinger boson representation. This nat- 
urally calls for a critical investigation of the SMA. To 
check the validity of SMA, we start with the mean 
field and find separate spin component mode functions 
(^a) = <i>Q (at zero temperature). The dynamics of $q, 
for the ground state is governed by Hs+Ha, which obeys 
the following coupled GP equation 




FIG. 2: The original (left column) and the renormalized 
(right column) wave functions along radial (upper panel) 
and axial (lower panel) directions for the -I- (solid line) and 
the — (dashed line) spin components. Other parameters are 
iV = 3.16 X 10^ X/iV = 0.5, and A = 2. All lengths are in 
units of y^h/mcOr- 

ih^o = [7i: + C2(n+ + ?!_)] $0 + 2c2<I>o*+'S-, (5) 
= [H-F6 + C2(n_ +no -n+)]$_ +C2$^$;, 

with n = -fi^V^/2M + Voxt + cqti, 12^ = and n = 

+ nQ + n_. We have developed a reliable numerical al- 
gorithm based on propagating Eq. (^) in imaginary time 
{it) that converges to the ground state while maintaining 
the conservation of both N and M . We take the initial 
wave function to be a complex Gaussian with a constant 



velocity, i.e., e' 



-(x'/2qi+y'/2q'+z'/2qi)-ik.f 



where Qx, Qy, 



Qz, and k are adjustable parameters that shall not affect 
the final converged solution. In the simplest case for the 
ground state, we assume $Q(r) = |<I'Q(r)|e'^° with 9a a 
global phase independent of r. Then only the relative 
phase A ~ 26q — 9^ — 9- shows up in with a term 
oc C2|<i>+$_$ol cos A. This gives A = (for C2 < 0) or tt 
(for C2 > 0) when is minimized p5[ , a conclusion also 
verified by numerical calculations. As first stated by Ho 
pT| , the spinor-1 condensate Hamiltonian H = Hs + Ha 
is invariant under gauge transformation e** and spin ro- 



tations U {a, P, t) 



iF^a -iFyfi -iF^T 



For the ground 



state that conserves however, the spin rotation sym- 
metry is reduced to the subgroup SO (2) generated by 
Thus irrespective of the signs of C2 , a transfor- 



mation of the form ^-^''o g-iFz{e+-6^) /2 ^^^^ always reduce 
a complex solutions to a real one pst . 

When i3 = as for ferromagnetic interactions with any 
values of magnetization < A'' or for anti-ferromagnetic 
interactions with = 0, we find = |(/>-| from the 
symmetry of Eq. (ph. We then rescale the wave func- 
tion (f)a — ^a/VNa such that (j)a is normalized to 1 
(J dr|$p(f)p = Nf^, the number of atoms in ^-th compo- 
nent), the asymmetric interaction energy then becomes 



(A+|<^^ 
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FIG. 3: M and dependence of the overlap integral 
|{0_|^SMA>| for ^^Na atoms (A = 1). |(0+|(?!>sma)| is essen- 
tially 1 to within ±0.001 in the same region. 



27V+7Vo|0+n0o|' + 2iV_iVo|0-n,/.o 
4iVo ( A^- ) 1 00 n 0+ 1 1 0- 1 cos A 



For ferromagnetic interactions (c2 < and A = 0), 
thus prove in general that Ea is minimized when 



10+1 = |0o| = 10-1 = 101, 



(6) 



we 



(7) 



and Nq/N = {\~M'^ /N'^)/2. The latter result (indepen- 
dent of all other parameters) was first obtained in Ref. 
|p3[ assuming the SMA, i.e. essentially assuming Eq. 

Our numerical solutions closely follow this as shown 
in Fig. |l|. For anti- ferromagnetic interactions (c2 > 0), 
S = holds only when = 0. In this case, using A = tt, 
we prove in general that is minimized to zero under 
Eq. (0), while No can be any value < 

For anti-ferromagnetic interactions (A^ ^ 0), we find 
that mode functions for the three spin components are 
different (see Fig. 2). Further analysis show that Ea is 
minimized if A'o = JT^ . 

We now discuss the relationship of Eq. (R) to the SMA 
Eq. (|^) . We note that the validity of Eq. (0) (including 
Ha) is in fact not equivalent to the validity of the SMA 
(excluding Ha)- For ferromagnetic interactions, with Eq. 
(0) and the relation between A^o and equation (||) 
simplifies to 



' 2M 



Kxt + (co+C2)iV|0|^ 



(f) = Ai0(f). (8) 



This shows that 0(r) is independent of M, and its de- 
viation from 0SMA comes only from the C2 term. This 
result can in fact be easily understood. Since cq + C2 = 
47r7i^a2/A/, 0(r) of Eq. (|^) is simply the ground state 
of the GP equation for an atomic scattering length of 
a2- In a ferromagnetic state, atomic spins are aligned 
locally. Two such atoms (Fi^2 = 1) only collide in 
the symmetric total spin F — 2 channel. For quanti- 
tative results, we compared |(0|0sma)| for ^''Rb atoms 
with ao — 101. Sob and 02 = 100. 4aB (a^ the 

Bohr radius). Other assumptions are: typical radial 
trap frequency = 27r x 10"^ (Hz), axial trap frequency 
Lu, = \LOr, and A = 0.1, 1, and 10. We also took 
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FIG. 4: The Ai dependence of energy components (in 
units of huir) (solid line), Ea (dashed line), and BM 
(dash-dotted line) for ^^Na atom at A = 1 and iV = 3.16 x 10^. 

N — 3.16 X 10^ and varied the total magnetization A4/N 
from to 1. Under these conditions, we find uniformly 
that |(0|0sma)| ~ 1 essentially because \c2\ \co\. 

For anti-ferromagnetic interactions (c2 > 0), two spe- 
cial cases arise: 1) when M = 0, using Eq. (j|), we prove 
that 1 0a I = |0sma|, which means the SMA is exact in 
this case; 2) when = A^, 0+ satisfies the same equa- 
tion as Eq. (||), so its derivation from the SMA only 
originates from the C2 term. For ^'^Na atoms, we use 
ao = 50 as and 02 — 55 as |2^ as an example in this 
case. Other parameters used are the same as in the fer- 
romagnetic case. Since iVo — 0, we consider only the 
± components. Figure 2 shows the original and renor- 
malized wave function for N — 3.16 x 10^, M/N — 0.5, 
and A = 2. We see clearly that (p^ and 0_ are different. 
Figure |3| shows the magnetization and atom number de- 
pendence of |(0_|0sma)| for a spherical trap. Since the 
-|- component contains the majority number of atoms, 
it is natural to find |(0+|0sma)| ~ 1- The value of 
|(0+|0sma)| at M N also indicates that the devia- 
tion contributed by C2 alone is also small for ^^Na atoms. 
While for |(0_|0sma)|, we see it becomes as low as 0.75 
when N = 3.16 x 10^ and when M approaches TV. To 
distinguish the different sources of deviations, we plot 
E^^ = (co/2) / drn"^, Ea, and BM in Fig. § We see 
that the BM term contributes the most. In Fig. ||, the 
overlap integral is shown to also depend on the trap as- 
pect ratio A. 

Finally, we discuss the implications of our result on the 
macroscopic alignment of the total spin of a spinor con- 
densate. For ferromagnetic interactions, the spatial dis- 
tribution of the total spin (F(f)) = J2ai3 ^H'^)^ai3^i3ir) 
is found to be pointed along the same direction, i.e. in- 
dependent of the spatial coordinates. Using A = 0, 
No = N{1 ~ M^/N^)/2, and iV± iV(l ± M/N)^/4:, 
it can be expressed as 



y/N^ - AJ2 cos( 



(F(r)) = |0(f)|^ ^VN-^ 



M'^ sin( 
M 



(9) 
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For anti- ferromagnetic interactions, we find 



(Fir)) = ( 

7V+|0+(f)|2-7V_|0_(f)|2 



(10) 



a state with all spins aligned in the ±z direction. It 
reduces to {F{r)) = for = 0. 
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FIG. 5: The overlap | ((/>±|^sma) | for A 
(dashed line), and 10 (dash-dotted line). 



0.8 



0.1 (solid line), 1 



To conclude, we presented both analytic and numeri- 
cal studies of the validity of the SMA. We find that de- 
viations of the ground state solution from the </)sma(^) 
come from two sources: the C2 or the B (due to conser- 
vation of M) term. For ferromagnetic interactions, the 



only source is the C2 term, which is negligible for ^^Rb 
atoms. One can therefore safely use the SMA. For anti- 
ferromagnetic interactions, if = 0, 0sma becomes the 
exact ground state wave function; For > 0, however, 
one can still use the (t>sMA for (/)_)-, but (f>- differs sig- 
nificantly if both N and Ai are large. In this case the 
BM term contributes the most to the deviation. Our 
conclusions from this study apply to the ground states 
of a spinor condensate. For dynamic problems Ref. [ p^ , 
the SMA may become worse. Our study suggests that 
instead of making the SMA as in Eq. (^), an improved 

ap</i^(r), where the mean- 



SMA could consist of = 



field solution $p(r) and its associated effective spin mode 
fimction 0^ = 

straints of conserved N and A4 . 



'N^i are obtained under the con- 



Such an approach can be 
important in studying beyond mean field quantum cor- 
relations. In a forthcoming article, we will report some 
results on condensate fragmentation. 

In summary, we have presented a detailed investiga- 
tion of the SMA for a spinor- 1 condensate and pointed 
out interesting structures of its ground state for both 
ferromagnetic and anti-ferromagnetic interactions. 

This work is supported by the NSF grant No. PHY- 
9722410 and by a grant from the National Security 
Agency (NSA), Advanced Research and Development 
Activity (ARDA), and the Defense Advanced Research 
Projects Agency (DARPA) under Army Research Office 
(ARO) Contract No. DAAD19-01- 1-0667. Partial sup- 
port from the NSF of China is also acknowledged. 



[1] See the extensive list of references 

[2: 



at; 



http: / / amo.phy.gasou.edu/bec.html/bibliography.htm^ 



[3] 



M. H. Anderson et aL, Science 269, 198 (1995); K. 
Davis et ai, Phys. Rev. Lett. 75, 3969 (1995); C. C. 
Bradley et al, ibid 75, 1687 (1995) and 79, 1170 (1997). 
Observation of Feshbach resonance has raised the hope 
that with a tunable asc, one can approach the interesting 



limit of nat 



1, even nai^ ^ 1 when exciting new 



[4] 

[5] 

[6] 

[7] 
[8] 

[9] 
[10] 

[11] 



physics may occur. See for example, Eric Braaten and 
H-W Hammer, Phys. Rev. Lett. 87, 160407 (2001). 
L.-M. Duan, J. I. Cirac, and P. ZoUer, |quant-ph/010705£ . 
H. Pu et al, Physica B 280, 27, (2000); C. K. Law et al., 
Phys. Rev. Lett. 81, 5257 (1998). 

E. Goldstein and P. Meystre, Phys. Rev. A 59, 3896 
(1998). 

T.-L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 
D. M. Stemper-Kurn et al, Phys. Rev. Lett. 80, 2027 
(1998). 

M. Barrett et al, Phys. Rev. Lett. 87, 010404 (2001). 
T. Ohmi and K. Machida, J. Phys. Soc. Jap. 67, 1822 
(1998). 

T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 



[12; 

[13 
[14 
[15 
[16 
[17 

[is: 

[19 

[20 
[21 
[22; 

[23; 

[24' 
[25 

[26 
[27 



C. V. Ciobanu et al, Phys. Rev. A 61, 033607 (2000). 
H. Pu et al, Phys. Rev. A 60, 1463 (1999). 
H. Pu et ai, Phys. Rev. A 61, 023602 (2000). 
H. Pu et al, Phys. Rev. A 63, 033603 (2001). 
S. Raghavan et al, Opt. Commun. 188, 149 (2001). 
M. Koashi and M. Ueda, Phys. Rev. Lett. 84, 1066 
(2000). 

T. Isoshima et ai, J 
M. Koashi and M. 

(2000) . 

M. Ueda, Phys. Rev. 
A. S. Sorensen et al, 
K. Helmerson and L 
(2001) 

U. V. Poulsen and K. Molmer, Phys. Rev. A 64, 013616 

(2001) . 

A. S. Sorensen, ( |cond-mat/0110372| ) 

T. Isoshima et ai. Phvs. Rev. A 6 0. 4857 (1999); N. P. 

Robins et ai, ( gond-matt /0105353 ) . 

E.G.M. van Kempen et ai, ( |cond-mat/01 106101 ). 

A. Crubellier et ai. Euro. Phys. Jour. D 6 (2), 211 (1999). 



. Phys. Soc. Jpn. 69, 3864 (2000). 
Ueda, Phys. Rev. Lett. 84, 1066 

A 63, 013601 (2000). 

Nature (London) 409, 63 (2001). 

You, Phys. Rev. Lett. 87, 170402 



